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Abstract 

The spiraling of adjacent trajectories in chaotic dynamical systems can be characterized by the 
distribution of local angular velocities of rotation of the displacement vector, which is governed 
by linearized equations of motion. This distribution, akin to that of local Lyapunov exponents, is 
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studied for three examples of three-dimensional flows. Toy model shows that the rotation rate of 
adjacent trajectories influences on the rate of mixing of dynamic variables and on the sensitivity 
of trajectories to perturbations. 
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I. INTRODUCTION 

In nonlinear dynamics of systems - flows with equations of motion 

x = F(x), (1) 

where x(t) is a vector representing the state of the system in a i^-dimensional phase space, 
x = {xi}, 1 < i < K, definitive role is played by the properties of evolution of the displace- 
ment vector r(t) that is governed by the linear system of equations 

r = Mr. (2) 

Here M is the local stability matrix with the elements M t j = dFi/dxj, that are taken at 
the points of the phase trajectory x(t) and depend, generally, on time. The existence of 
the exponentially growing solutions of the system (2) serves as a definition of chaoticity of 
motion of the system (1). 

For almost any initial conditions the displacement vector r(t) will approach the eigenmode 
of the system (2) that has the largest rate of the exponential growth, given by the first 
(maximal) characteristic Lyapunov exponent 0\ or simply the Lyapunov exponent a [l|, |2| . 
Thus the behaviour of r(t) at large times t ^> a -1 loses its dependence on the initial 
conditions and represents some characteristics of the chaotic component as a whole. This 
property serves as a base for the standard method of numerical calculation of the Lyapunov 
exponent [3| , that is determined from the rate of growth of the length of vector r (t) . However, 
the evolution of the orientation of this vector is usually ignored. 

We restrict ourselves by the phase space with dimensionality K = 3, the minimal one in 
which the chaotic motion of systems - flows like (1) is possible. Then the evolution of the 
direction of the displacement vector r(t) can be described as a rotation around the direction 
of the phase trajectory: the adjacent trajectories can be spiraling, as it is shown in figure 1. 




FIG. 1: Qualitative scheme of position of adjacent phase trajectories for chaotic motion in three- 
dimensional phase space. Whereas the displacement vector grows in length by the exponential law 
with the rate a, it may also rotate with the angular velocity oo 



The instantaneous angular velocity of this rotation uj{t) (we shall also call it "the rotation 
rate") is denned by the expression 

1 r(£ + e) x r(t) • v(t) 



u(t) 



lim 



(3) 



|r(* + e)|| 

where v(t) = x(£) is a vector of phase velocity, that is tangent to the phase trajectory, and 
a x b • c denotes the mixed product of three vectors. 

The distribution W of values of u(t) does not depend on the choice of the initial moment 
of time or the initial conditions for the system (2) as far as they belong to the same chaotic 
component. It is as universal as the distribution of the local Lyapunov exponents that are 
defined by the expression 

a(t) = limilnM±£M- (4) 



The latter quantity is widely used for studies of the structure of the phase space of chaotic 
systems |4|, |s|, |6(. From comparison of equations (J3J) and (jlj) one can assume that the 
instantaneous angular velocity in Pickwickian sense can be considered as an imaginary part 
of the local Lyapunov exponent and could be used as a probe for exploration of fine details 
of chaotic dynamics. 



II. THREE EXAMPLES 



The distribution W{u) could be easily obtained numerically. As the first example we 
consider the famous Lorenz model jj| given by the equations of motion 



X = -aX + aY, 

Y = -XZ + rX-Y, 

Z = XY - bZ, 



(5) 



with the standard values of parameters a = 10, r = 28, b = 8/3. The distribution W{ui) is 
shown in figure 2. 
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FIG. 2: Distribution W of the instantaneous rotation rate ui for the Lorenz model defined by the 
system of equations (5) with the standard values of parameters a = 10, r = 28, b = 8/3. 



The average value of the rate of rotation u = 0.158 is much less than the width of the 
distribution (standard deviation of a;) Au; = 0.63. It may be noted that for the Lorenz model 
with the used values of parameters the characteristic frequency of motion Q estimated from 
the power spectrum of the variable Z(t) is about Q ~ 8, whereas the Lyapunov exponent 
a w 1.0. Thus the small value of u may indicate the presence of the additional time scale 
of motion. We also note that the time of decay of correlations of u(t) is rather close to 



T 



For the second example we take a conservative autonomous system, the Pullen - Edmonds 



model oscillator 81 with the Hamiltonian 



H = y + U(x, y)= l - (pi +P 2 y )+\ (x 2 + y 2 + x 2 y 2 ) , (6) 

where x,y are the Cartesian coordinates of the particle and p x ,P y are the corresponding 
canonically conjugated momenta. In expression (6) the particle mass m, the frequency 
of small oscillations ujq and the nonlinearity length A have been used as the unit scales. 
The motion of this system takes place on a three-dimensional energy surface in the four- 
dimensional phase space. Since the energy of the system is conserved, E = H = const, the 



absolute value of the total momentum p (x, y) = w 2 [E — U (x, y)} is completely defined by 
the position of the phase point in the configurational space. The angle tp that gives the 
direction of the total momentum is defined by the equations 



p x =p cos ip, p y =p sin if. 
It can be taken as the third dynamic variable, that yields the equations of motion 



1 (dU 



dU 



x = pcos(p, y=psm.(p, (p = - — — sin (p — — COS(y9 

p \ ox ay 



(7) 



(8) 



For the motion on the energy surface E = 15, that is nearly ergodic J2|, the distribution 
W(oS) is shown in figure 3. 
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FIG. 3: Distribution W of the instantaneous rotation rate uj for the Pullen - Edmonds model, 
defined by the system of equations (8), on the energy surface E — 15. 



The average value of the rate of rotation uj = —1.49 is comparable to the width of the 
distribution Auj = 1.53 and to the characteristic frequency of motion Q ~ 1. 



For the third example we take a conservative non-autonomous system with the equation 
of motion 



x + x 



F cos Qt, 



(9) 



that describes an oscillator with cubic nonlinearity under the influence of external har- 
monic force (Duffing model). The system (9) could be turned into an autonomous three- 
dimensional system by introducing dynamical variables x, y = x and interpreting the time 
t in the RHS as a third dynamic variable z that obeys the equation of motion z = 1 |l( . For 
the motion in the chaotic component that surrounds the line x = 0, y = the distribution 
W(oS) is shown in figure 4. 
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FIG. 4: Distribution VF of the instantaneous rotation rate u for the Duffing model defined by 
equation (9) with parameters F = 5, O = 1. 

It follows from equation (9) that the rotation rate is always negative. Its average value 
oo = —2.24 is comparable to the width of the distribution Auj = 2.38. We note that the 
value of io is not trivial: it is incommensurate with the driving frequency Q = 1. 

From comparison of graphs in figures 2 - 4 it may be noticed that although our examples 
were taken from very different classes of dynamical systems, the forms of distributions have 
some common features: they are asymmetric, strongly peaked at a single non-zero most 
probable value and have extended wings. 



III. SPIRALING AND MIXING 

The proof of usefulness of the quantity uj(t) for the exploration of the phase space we 
leave for the future studies. There is another way to reveal the relevance of u(t) to the 
properties of chaotic motion: to find a relation between this quantity and some established 
characteristics of chaoticity. The affinity between uo and a suggests that relations that 
include the Lyapunov exponent may be extended to include u in a similar way. 

In nonlinear dynamics there is a widely known rule of thumb stating that the rate of 
mixing (the inverse time of damping of correlations of dynamic variables) 7 is approximately 
equal to the Lyapunov exponent: 

7Riff. (10) 

This relation sometimes turns into exact equality (e.g. for the linear random number gener- 
ator x' = {Qx}, where {, } denotes the fractional part of a number, with integer parameter 
Q); sometimes equation (10) reflects the scaling property 7 oc a (e.g. it holds exactly for 
billiards), and in general it is expected to be fulfilled at least semiquantitatively for arbi- 
trary systems [HJ] . From the picture in figure 1 one could surmise that the increase of u will 
intensify the rate of exploration of the phase space and thus enhance the mixing. 

For the check of this hypothesis we retreat to two-dimensional mappings, that can be 
interpreted as Poincare mappings for three-dimensional flows. In this class it is easy to 
construct a model that permits controllable and independent changes of to and a. Let 
{£ n , r) n } be the components of the vector s n of normalized displacement at the n-th iteration, 
s n = r n /||r n ||. Then the rate of rotation is equal to the angle of rotation for one step: 

uj n+1 = arcsin (£ n +ir)n - £ n r)n+i) ■ (H) 

Let's define the rate of damping of correlations (rate of mixing) of a dynamical variable 
with the known correlation function B(n) by the relation 

^(s3l)' s -s(o)gl B <">l < 12 > 

For the exponentially decreasing correlation function B (n) = B (0) exp (—an) this definition 
produces the rate of mixing equal to the exponent, 7 = 0. 

For the study of influence of rotation on mixing we take a two-dimensional mapping for 
which it is possible to variate rotation rate u with constant a, namely a piecewise linear 



mapping of a unit square on itself 



x' = {ax + by} , y' = {—bx + ay} , 



(13) 



where 

a = expacosu, b = expo" sin u. (14) 

The behaviour of rates of mixing for dynamical variables x and y is shown in figure 5. 
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FIG. 5: Dependence of the mixing rate 7 on the rotation rate u> for the two-dimensional mapping 
(13): (a) - for variable x; (6) - for variable y. 

The results of numerical calculation confirms the intuitive guess: at small values of uo its 
increase leads to the increase of the rate of mixing 7. 



IV. SPIRALING AND RESPONSE TO PERTURBATIONS 



Let's compare the system (1) to its perturbed version with the equations of motion 



x = F f x) + eV (x) 



(15) 



where eV(x) is some static perturbation. Let xo(t) and xy(t) be the laws of motion of the 
unperturbed system (1) and perturbed system (15) correspondingly with the same initial 
conditions, and let A (t) = |x (t) — xy (t)\ denote the distance between the phase trajecto- 
ries at a given moment t. For a given point of the phase space let's define a quantity 

A(x)= lim ^expHxt), (16) 

that defines the amplitude of the response of the system to the perturbation of a given form. 
This quantity, averaged over a given chaotic component of the phase space, A = (A (x)), 
gives a convenient measure of sensitivity of the system (1) to small perturbations. 

Let's treat the simplest model of the evolution of deviations given by a system of flows 

i - <r£ + ujrj = U, i] - ui - arj = V, (17) 

with constant parameters a, u, U and V . In the absence of perturbations (U — V — 0) it 
describes the rotation of the displacement vector with the constant angular velocity uo and 
its exponential growth with the exponent a. For U, V ^ the asymptotics of solutions for 
er£ 3> 1 have the form 

e at 

£ ss — [(all - ujV) cos ut + (uU + aV) sin out] , 

a 2 + u 2 

at 

r] sa — \(uU + aV) cos cut + (-all + uV) sin cot] . (18) 

a 2 + uj 2 

By definition (16) for each of the perturbations we have 

A= l (19) 

<2 (a 2 + u 2 ) 



Therefore the schematic model (17) shows that with the increase of the rotation rate the 
amplitude of the response will decrease. One may infer that this connection may be universal. 
Let's return to the piecewise model (13) perturbed by small variations of its control 
parameters: 

x = {{a + e)x + by}, y' = {-bx + ay} (20) 

and 

x' = {ax + (b + e) y} , y' — {—bx + ay} . (21) 

Fig. 6 demonstrates that with the increase of the rotation rate u the averaged amplitude 
of response to the perturbation A in general decreases, albeit non-monotonously, in analogy 
with the dependence (19). 
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FIG. 6: Dependence of the response amplitude A on the rotation rate u for different perturbations 
of the two-dimensional mapping (13): (a) - perturbed system is given by equation (20); (b) - 
perturbed system is given by equation (21). In both pictures filled circles are for a = 0.25 and 
open circles - for a = 0.5. 

Recently the changes of the evolution under the influence of small perturbation have 
attracted considerable attention in the theory of Hamiltonian chaotic systems in connection 
with the problem of dynamics of fidelity (llj, [l2j, [l3j, 1 141 ]. The example described in this 
section shows that the rotation rate can display itself in the preexponential factors of this 
response. 



V. GENERALIZATION TO HIGHER DIMENSIONS 



Everywhere above we have restricted ourselves to analysis of trajectories in three- 
dimensional phase space. The generalization of the definition (3) to higher dimensions 
could be obtained straightforwardly from the expression for the oriented volume Vk in K- 
dimensional space [l^J. It can be constructed from the unit vector along the direction of 
phase velocity u(t) = v(t)/||v(t)|| and (K — 1) unit vectors along the displacement vec- 
tors taken in consequent equidistant moments of time, s(t + ke) = r(£ + A;e)/||r(t + ke)\\, 
1 < k < K — 2 by their convolution with the completely antisymmetric tensor of the i^-th 
rank E ijk ,„ n : 



V K (t) = Enk... n UiSj (t) s k (t + e) ....s n {t + {K- 2) e) . (22) 

However, this quantity scales as e N , where N = (K — 1)(K — 2)/2, and the generalized rate 
of "rotation", 

u K {t) = \ime- N V K {t), (23) 

for K > 3 depends not only on the elements of stability matrix at a given moment of time 
but also on their time derivatives up to the order N — 1. 
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